Vortex states in superconductors with strong Pauli-paramagnetic effect 
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Using the quasiclassical theory, we analyze the vortex structure of strong-paramagnetic super- 
conductors. There, induced paramagnetic moments are accumulated exclusively around the vortex 
core. We quantitatively evaluate the significant paramagnetic effect in the .ff-dependence of vari- 
ous quantities, such as low temperature specific heat, Knight shift, magnetization and the flux line 
lattice (FLL) form factor. The anomalous H-dependence of the FLL form factor observed by the 
small angle neutron scattering in CeCoIns is attributable to the large paramagnetic contribution. 

PACS numbers: 74.25.Op, 74.25.Ha, 74.25.Jb, 74.70.Tx 



I. INTRODUCTION 

Recent two seemingly quite independent activities 
prompt us to cope those in a unified way because two 
typical experiments in each field suggest the vortex struc- 
ture associated with strong influence of the mismatched 
Fermi surface, namely two Fermi levels for up and down 
spins are different; One is in the rotating fermion su- 
perfluids of neutral 6 Li atom gases under the popu- 
lation imbalanceji^^ii. The other is heavy fermion su- 
perconductors with Zeeman-shifted Fermi surfaces un- 
der an applied field H. A heavy fermion compound 
CeCoIns is a prime candidate of a superconductor with 
strong Pauli-paramagnetic effect, where at higher fields 
the upper critical field H C 2 changes to the first or- 
der phase transitional 7 - and new Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state appears4i2.I2.iiiI2.i2.IliS, In 
the FFLO states, the pair potential is considered to have 
periodic spatial modulation in addition to the vortex 
modulation, 16 . 17 d8il9 1 20,2i 1 22 1 23 1 24,25 It is be cause Cooper 

pairs of up- and down-spins acquire non-zero momentum 
for the center of mass coordinate of the Cooper pair by 
the Fermi surface splitting of up- and down-spin electron 
bands due to Zeeman shift. The Fermi surface splitting 
is also an origin of the Pauli-paramagnetic pair breaking 
and the appearance of paramagnetic moments. There- 
fore, even when the vortex states do not enter to the 
FFLO state yet, the strong paramagnetic effects may se- 
riously contribute to the vortex state in superconductors. 

It is expected that, in the presence of strong param- 
agnetic effect, observed quantized vortices in both sys- 
tems of atomic gases and solid states should have univer- 
sal common properties that are absent in conventional 
vortex picture. As for the rotating atomic gas under 
population imbalance within a trap potential, the vortex 
states were studied by Bogoliubov-de Gennes theory in 
the configuration of single vortex4. There, paramagnetic 
moments are enhanced in the vortex core region, and 
the vortex core structure is related to the spectral evolu- 
tion of quasiparticles around the vortex in the presence 
of Zeeman shift. Thus, also in solid states it is neces- 
sary to study the vortex states in superconductors with 
strong paramagnetic effect, clarifying the exotic vortex 
core structure of pair potential, paramagnetic moments, 



internal magnetic field distributions, and local electronic 
states. These paramagnetic effects also give significant 
contributions to bulk properties, such as specific heat, 
paramagnetic susceptibility, or magnetizations. 

In some heavy fermion superconductors, the param- 
agnetic effects due to Zeeman shift are important to 
understand the properties of the vortex states, because 
the superconductivity survives until under high magnetic 
fields due to the effective mass enhancement. The para- 
magnetic contributions are eminent at higher fields. For 
example, the .ff-dependence of low temperature specific 
heat C(H) is often used to distinguish the presence of 
nodes in the pairing potential. The Sommerfeld coeffi- 
cient j(H) — liniT^o C(H)/T cx H in s-wave full gap 
superconductors, and j(H) cx y/H by the Volovik effect 
in d-wave pairing with line nodesi 26 : 27 ' 28 : 29 The curves 
of j(H) are expected to smoothly recover to the normal 
state value towards H C 2. However, in some heavy fermion 
superconductors, C{H) deviates from these curves. In 
CeCoIn 5 , C{H) shows concave curves, i.e., C{H) oc H a 
(a > I) at higher fields.— i2L22 This behavior is difficult 
to be understood only by effects of the pairing symme- 
try. A similar C(H) behavior is observed also in UBei3j22 
The experimental data of magnetization curve Mt ota \{H ) 
in CeCoIns show a concave curve at higher fields, in- 
stead of a conventional convex curve* 7 - As an anomalous 
behavior of CeCoIns, the small angle neutron scattering 
(SANS) experiment reported the ii-dependence of flux 
line lattice (FLL) form factor determined from the Bragg 
intensity. 34 While the form factor shows exponential de- 
cay as a function of H in many superconductors, it keeps 
almost constant at lower fields for H || c in CeCoIns. 
As for properties of CeCoIns, the contribution of anti- 
ferromagnetic fluctuation is also proposed in addition to 
the strong paramagnetic effect i 15 : 35 ' 36 : 37 Therefore, it is 
expected to study whether properties of vortex states in 
CeCoIn 5 are theoretically explained only by the param- 
agnetic effect. Theoretical studies of the 77-dependences 
also help us to estimate strength of the paramagnetic 
effect from experimental data of the if-dependences in 
various superconductors. 

In this paper, based on the selfconsistent micro- 
scopic calculation of quasiclassical Eilenberger the- 
ory , 27 i 28 i 29 i 38 i 39 i 40 we study the spatial structure of the 
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vortex states in the presence of the paramagnetic ef- 
fect by Zeeman-shiftj 25 i 41 ' 42 i 43 We will clarify the para- 
magnetic effect on the vortex core structure, calculating 
the pair potential, paramagnetic moment, internal mag- 
netic field, and local electronic states. We also study 
the paramagnetic effect by quantitatively estimating the 
-ff-dependence of low temperature specific heat, Knight 
shift, magnetization and FLL form factors, and show the 
relation of the -ff-dependence behaviors and the strength 
of paramagnetic effect. The anomalous field-dependence 
of FLL form factor observed by SANS experimental is 
explained by the strong paramagnetic effect. 

Previous works in the selfconsistent quasiclassical cal- 
culation were mainly applied to the vortex state in the 
absence of paramagnetic effects, and successfully esti- 
mate local electronic states around the vortex core and 
the ff-dependence of the low temperature specific heat, 
studying the effect of the pairing symmetry ) 27 i 28 ' 44 gap 
anisotropy^ and multi-band structured The selfcon- 
sistent calculation of the pair potential is necessary for 
quantitatively valid evaluation of the vortex states, since 
we have to use the correct vortex core size in the cal- 
culation at each field and temperature. As for previous 
quasiclassical studies on the strong paramagnetic effect 
on the vortex state, there were some works focusing on 
the FFLO vortex states^ii 2 . Without FFLO modulation 
but with strong paramagnetic effect, many studies were 
done along H C 2(T), and there were few studies far from 
H C 2. The quasiclassical study on the -ff-dependence of 
the specific heat and magnetization were done by Adachi 
et al. at T = 0.4T C , using Landau level expansion^ 1 - In 
this paper, we study the -ff-dependence at T = 0.1T C , 
solving the Eilenberger equation numerically by the ex- 
plosion metho d 39 ' 40 in the vortex lattice state, and also 
study the vortex core structure including local electronic 
states, and the -ff-dependence of paramagnetic suscepti- 
bility and the FLL form factors. 

After giving our formulation of quasiclassical theory 
in Sec. HH we study the paramagnetic effect on the ff- 
dependence of paramagnetic susceptibility and low tem- 
perature specific heat in Sec. IIIIl In Sec. IIV1 we show the 
paramagnetic contributions on the vortex core structure, 
and the local electronic state in the presence of Zeeman 
shift. In Sec. [V] we estimate the -ff-dependence of FLL 
form factor, and discuss the anomalous .ff-dependence 
observed by SANS in CeCoIns. The paramagnetic effect 
on the magnetization curve is presented in Sec. IVI1 and 
the last section is devoted to summary and discussions. 



II. QUASICLASSICAL THEORY INCLUDING 
PARAMAGNETIC CONTRIBUTION 

We calculate the spatial structure of the vortex lat- 
tice state by quasiclassical Eilenberger theory in the 
clean limit ) 27 i 28 ' 29 ' 38 ' 39 i 40 including the paramagnetic ef- 
fects due to the Zeeman term /-tB-B(r), where B(r) is 
the flux density of the internal field and fiB is a renor- 



malized Bohr magneton i 25 ' ' ' The quasiclassical the- 
ory is quantitatively valid when £ 3> 1/fcp (&f is the 
Fermi wave number, and £ is the superconducting coher- 
ence length), which is satisfied in most of superconduc- 
tors in solid states. The quasiclassical Green's functions 
g(u n +ifj,B,k, r), f(w n +ifi,B, k, r), and f f (uj n +ifiB, k, r) 
are calculated in the vortex lattice state by the Eilen- 
berger equation 

{u n + i/j,B + v • (V + iA)}/ = Afa, 

{oj n + i/i-B - v • (V - iA)} p = A*cj ) *g, (1) 

where g — (1 — /J^) 1 ^ 2 , Re<? > 0, v = v/ufo, and 
fi = fiBBo/irkBT c . We mainly consider the (i-wave pair- 
ing case for a pairing function, </>(k) = \/2cos2#, as sug- 
gested in CeCoIns^i 3 - 1 - k is the relative momentum of the 
Cooper pair, and r is the center-of-mass coordinate of the 
pair. We consider the case of two-dimensional cylindrical 
Fermi surface, k = kp(cos9,sm9), where < 9 < 2ir. 
The Fermi velocity is given by v = t>po(cos 9, sin#). 
Throughout this paper, length, temperature, and mag- 
netic field are scaled by Rq, T c , and Bq, respectively. 
Here, i?o = hvF /2irkBT c , B = hc/2\e\B,Q. Matsubara 
frequency uj n = (2n + l)irT, energy E, and pair potential 
A are scaled by irkBT c . Since magnetic fields are applied 
to the z axis direction, in the symmetric gauge the vec- 
tor potential A(r) = x r + a(r), where B = (0, 0, B) 
is a uniform flux density and a(r) is related to the in- 
ternal field B(r) = B + V x a(r). The unit cell of 
the vortex lattice is given by r = si(ui — 112) + S2U2 
with -0.5 < s l < 0.5 (i=l, 2), Ui = (a, 0,0) and 
u 2 = (o/2, a y , 0). In the <i-wave pairing, we consider the 
case of square vortex lattice, i.e., a y /a = 1/2, where the 
nearest neighbor vortices are located to the nodal direc- 
tion. In the d-wave superconductors, this square lattice 
configuration is stable at higher fields i 27 ' 34 

The pair potential is selfconsistently calculated by 

A(r) =g Q N Q T £ (^(k) (/ + /t*)) k (2) 

0<w„<w cut 

with (.goiVo)- 1 = lnT + 2TE < Wn < Wcut w" 1 . (■ ■ -) k indi- 
cates the Fermi surface average. We use w cu t = 20/cbT c . 
The vector potential for the internal magnetic field is 
selfconsistently determined by 

2T 

V x (V x A) = V x M para (r) - — V (v F Im ff ) k , (3) 

0<UJ n 

where we consider both the diamagnetic contribution of 
supercurrent in the last term and the contribution of the 
paramagnetic moment M para (r) = (0, 0, M para (r)) with 

Af para (r) = M (^pL E < Im M>k) ■ (4) 

The normal state paramagnetic moment Mq = {p/k) 2 B, 
k = B Q /nkBT c y/8TTN and N is the density of states 
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(DOS) at the Fermi energy in the normal state. We 
solve Eq. |T]) and Eqs. ©-© alternately, and obtain 
selfconsistent solutions as in previous works , 27 i 40 i 45 un- 
der a given unit cell of the vortex lattice. The unit cell is 
divided to 41 x 41 mesh-points, where we obtain the qua- 
siclassical Green's functions, A(r), M para (r) and A(r). 
When we solve Eq. |T]) by the explosion method, we 
estimate A(r) and A(r) at arbitrary positions by the in- 
terpolation from their values at the mesh points, and by 
the periodic boundary condition of the unit cell including 
the phase factor due to the magnetic field . 27 i 39 i 40 i 45 

Using Doria-Gubernatis-Rainer scaling ) 43 ! 46 we obtain 
the relation of B and the external field H as 



ot£. 2 \ g + l 

+uiRc{g - l}) k )r, 



(5) 



where (• • -) r indicates the spatial average. We con- 
sider the case of large Ginzburg-Landau (GL) parameter 
kgl ~ k = 89 and low temperature T/T c = 0.1. For the 
two-dimensional Fermi surface, k = (7C(3)/8) 1 / 2 kql ~ 
kgl^ 7 - In these parameters, \B — H\ < 10~ 4 £>o- 

When we calculate the electronic states, we solve Eq. 
(fT]) with \u> n — ► E+irj. The local density of states (LDOS) 
is given by N(r, E) = 7V T (r, E) + Ni(r, E), where 

N a (r,E) = N (Rc{g(Lu n + icTnB,k,r)\ lL , n ^E+ i r 1 })k (6) 

with a — 1 (—1) for up (down) spin component. We 
typically use r\ = 0.01. The DOS is obtained by the spa- 
tial average of the LDOS as N(E) = N^E) + N^E) = 
(N(r,E)) r . 



III. FIELD DEPENDENCE OF 
PARAMAGNETIC SUSCEPTIBILITY AND 
ZERO-ENERGY DOS 

First, we discuss the field dependence of zero-energy 
DOS y(H) — N(E = 0)/No and paramagnetic suscep- 
tibility x{H) = (A/ para (r)) r /Afo, which are normalized 
by the normal state values. From low temperature spe- 
cific heats C, we obtain ~f(H) cx C/T experimentally. 
And x{H) is observed by the Knight shift in NMR ex- 
periments, which measure the paramagnetic component 
induced by an external field via the hyperfine coupling 
between a nuclear spin and conduction electrons. In d- 
wave pairing, x{H) shows v^ff-behavior at low fields4^ 

As shown in Fig. [TJ 7 (dashed lines) and % (solid 
lines) show almost the same behavior at low temper- 
atures. First, we see the case of d-wave pairing with 
line nodes in Fig. [Ha). There j(H) and x(-ff) describe 
\f~H-\\ks recovery smoothly to the normal state value 
(7 — X = 1 & t ^2) m the cas e of negligible param- 
agnetic effect (n = 0.02)j 26 i 27 i 28 i 29 With increasing the 




FIG. 1: (color online) The magnetic field dependence of the 
paramagnetic susceptibility x(^) (solid lines) and the zero- 
energy DOS j(H) (dashed lines) at T = 0.1T C for various 
paramagnetic parameters fi = 0.02, 0.86, 1.7, and 2.6 in the 
d-wave (a) and s-wave (b) pairing cases. The insets show 
same data as a function of H/H C 2- 



paramagnetic parameter /i, H C 2 is suppressed and the 
Volovik curve 'j(H) oc gradually changes into curves 
with a concave curvature. For large jj,, H C 2 changes to 
first order phase transition^ We note that at lower fields 
all curves exhibit a \fli behavior because the paramag- 
netic effect (oc H) is not effective. Further increasing H, 
7 (J?) behaves quite differently. There we find a turning 
point field which separates a convex curve at lower H 
and a concave curve at higher H. This inflection point 
increases as /i decreases. In the inset of Fig. [T] we plot 
'y(H) and x{H) as a function of H/H C 2- The overall 
-H-dependence at < H < H C 2 can be used to analyze 
experimental data, in order to estimate the strength of 
the paramagnetic effect, /1. 

To examine effects of the pairing symmetry, we show 
j{H) and x(H) a ls° f° r s-wave pairing 0(k) = 1 in 
Fig. [ljb), where we use the triangular lattice configu- 
ration a,y/a = V3/2. In the independence of j(H) and 
x(H), differences by the vortex lattice configuration of 
triangular or square are small. The difference in the in- 
dependences of Figs. [Ua) and [TJb) at low fields comes 
from the gap structure of the pairing function. In the 
full gap case of s-wave pairing, j(H) and x(H) show 
ii-linear-like behavior at low fields. With increasing the 
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FIG. 2: (color online) Spatial structure of the pair potential 
(a), paramagnetic moment (b) and internal magnetic field 
(c) at T = O.lTc and B = O.LB , where a = 11.2i? , in 
d-wave pairing. The left panels show |A(r)|, M para (r), and 
B(r) within a unit cell of the square vortex lattice at /i = 1.7. 
The right panels show the profiles along the trajectory r from 
the vortex center to the midpoints between nearest neighbor 
vortices, fx = 0.02, 0.86, 1.7, and 2.6. 



paramagnetic effect, iJ-linear behaviors gradually change 
into curves with a concave curvature. As seen in Figs. 
[IJa) and[ljb), paramagnetic effects appear similarly at 
high fields both for s-wave and d-wave pairings. 

The concave curves of the specific heat at higher fields 
by strong paramagnetic effect were also presented in Ref. 
l4ll at T = 0.4T C . In the present calculation, the concave 
curves are confirmed even at low temperature T = 0.1T C , 
where C/T ~ j(H) without large specific heat jump at 
H C 2. In j(H) at low T, the differences of s-wave and 
d-wave pairings at lower fields are clearly seen. 



IV. PARAMAGNETIC CONTRIBUTION ON 
VORTEX CORE STRUCTURE 

In order to understand the contribution of the param- 
agnetic effect on the vortex structure, we illustrate the 
local structures of the pair potential |A(r)|, paramag- 
netic moment M para (r), and internal magnetic field B(r) 
within a unit cell of the vortex lattice in Fig. O Since we 
assume d-wave pairing with the line node gap here, the 
vortex core structure is deformed to fourfold symmetric 
shape around a vortex core. 27 ' 49 It is noted that the para- 
magnetic moment is enhanced exclusively around the 
vortex core, as shown in Fig. [2Jb) where the four ridges 



FIG. 3: (color online) Local density of states at r/Ro = 
(a), 0.8 (b) and 1.6 (c) from the vortex center towards the 
nearest neighbor vortex direction in d-wave pairing. Solid 
lines show JVf(r, E)/No for up-spin electrons, and dashed lines 
show N L (r,E)/N at B = 0.1B . fj, = 1.7 and T = 0.1T C . 
(d) Spatial-averaged DOS at B/Bo = 0.1 and 0.01 in d-wave 
pairing. Solid lines show JVj (E) /No for up-spin electrons, and 
dashed lines show Ni(E)/No- 

of paramagnetic moment are extended towards the anti- 
nodal directions from the core. Since the contribution of 
the paramagnetic vortex core is enhanced with increas- 
ing /i, internal field B{v) consisting of diamagnetic and 
paramagnetic contributions is further enhanced around 
the vortex core by the paramagnetic effect, as shown in 
Fig. HJc). When [i is large, the pair potential |A(r)| is 
slightly suppressed around the paramagnetic vortex core, 
and the vortex core radius is enlarged, as shown in Fig. 
Ha). 

The enhancement of M para (r) around vortex core is 
related to spatial structure of the LDOS N a (r,E). As 
shown in Fig. [3ja), the LDOS spectrum shows zero- 
energy peak at the vortex center, but the spectrum is 
shifted to E — zk/iB due to Zeeman shift. The peak 
states at E > is empty for N^(E,r), and the peak 
at E < is occupied for Ni(E,r). Therefore, from the 
relation 

Mpara(r) = -/i B / (Aff (E, r) - N l (E,r))dE, (7) 

J — oo 

large M para (r) appears due to the local imbalance of up- 
and down-spin occupation around the vortex core. As 
shown in Figs. [3Jb) and[3£c), increasing the distance 
from the vortex center, the peak of the spectrum is split 
into two peaks, which are shifted to higher and lower 
energies, respectively. When one of split peaks crosses 
E = 0, the imbalance of up- and down-spin occupation is 
decreased. Thus, M para (r) is suppressed outside of vortex 
cores. 

In Fig. [3Jd), we present the spectrum of spatially- 
averaged DOS. In the DOS spectrum, peaks of the LDOS 
are smeared by the spatial average. Because of the flat 
spectrum at low energies, paramagnetic susceptibility 
x(H) shows almost the same iJ-behavior as the zero- 
energy DOS 7(i?) ~ N{E — 0) even for large fi, as shown 
in Fig. [TJ while x{H) counts the DOS contribution in the 
energy range \E\ < fxH, i.e., X (H) ~ tf" N t (E)dE/»H 
from Eq. 0. 
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FIG. 4: (color online) Field dependence of FLL form factor 
Fi, for ^ = 0.02, 0.86, 1.7, and 2.6 at T = 0.1T C in d-wave 
pairing, (a) |-Fi,o| 2 is plotted as a function of H. The ver- 
tical axis is in logarithmic scale, (b) We plot ln|Ji,o| 2 as a 
function of H/H C 2, where Fi,o is scaled so that H C 2 for each 
fj, corresponds to 5 [T] . Open circles are experimental data in 
CeCoIn 5 observed by SANS experiment.— 



V. FIELD DEPENDENCE OF FLUX LINE 
LATTICE FORM FACTOR 

One of the best ways to directly see the accumulation 
of the paramagnetic moment around the vortex core is 
to observe the Bragg scattering intensity of the FLL via 
SANS experiment. The intensity of the (h, fc)-diffraction 
peak is given by I h ,k = |-Ph.,fc| 2 /|qh.,fc| with the wave 
vector q/,.^ = h<\i + fcq 2 , qi = (27r/a, — 7r/a y , 0) and 
q2 = (2n/a,Tr/a y ,0). The Fourier component F^.k is 
given by B(r) = Y,h,k F h,kexp(iq htk ■ r). The intensity 
of the main peak at (h, k) = (1, 0) in the SANS for FLL 
observation probes the magnetic field contrast between 
the vortex cores and the surrounding. 

We calculate the field dependence of l-F^o | 2 , which is 
shown in Fig. 2J In the case of negligible paramagnetic 
effect (/i = 0.02), |-Fi,o| 2 decreases exponentially as a 
function of H . This is a result for the d-wave pairing at 
low T in the clean limit. With increasing paramagnetic 
effect, the decreasing slope of |-F\,o| 2 becomes gradual, 
and changes to increasing functions of H at lower fields 
in extremely strong paramagnetic case (/j, = 2.6). This is 
because |-Pi,o| includes enhanced paramagnetic contribu- 
tion proportional to fiH , reflecting the enhanced internal 



field around the vortex core, shown in Fig. [UJc), by the 
paramagnetic moment. 

The SANS experiment in CeCoIns reported that 
|-Fi,o| 2 is almost constant as a function of H within the 
field range 0.08 < H/H c2 < 0.4 (0.4[T] < H < 2.0[T]) 
for H || c. 34 This behavior is reproduced by our calcu- 
lation for /i ~ 1.7. There, |-Fi,o| 2 shows flat behavior at 
low fields, since the paramagnetic contribution increasing 
with H compensates the conventional decrease of |.Fi,o| 2 
as a function of H. For the comparison to the exper- 
imental data, we plot ln|Fi.o| 2 as a function of H/H C 2 
at lower fields in Fig. HJb). There, unit of the magnetic 
field in the calculated data for each fj, is rescaled so that 
H c2 corresponds to 5 [T], i.e., H c2 in CeCoIn 5 for H \\ c. 
For quantitative accordance of the results for /i = 1.7 
with the experimental data, we tune the GL parameter 
as h = 89. The variations of internal fields are roughly 
proportional to k -2 , as seen from Eq. Q- Changing 
k, we can shift curves in Fig. Ufb) towards the vertical 
direction. The slopes of the curves in Fig. H£b) are de- 
termined by the paramagnetic effect. With increasing /i, 
the negative slope becomes gradual, and even changes to 
positive slope. When the paramagnetic effect is negligi- 
ble (/i = 0.02), In |Fi,o| 2 decreases by 2 in the field range 
0.1 < H/H C 2 < 0.5, which corresponds to the exponen- 
tial decay as a function of H, as expected in conventional 
sup erconductors (also see the theoretical curves in Ref. 
[341 ). For the large paramagnetic case fx — 1.7, ln|i*i i0 | 2 
does not decrease as a function of H, which accords with 
the experimental data [circles in Fig. 0Jb)] of SANS ex- 
periments. The anomalous if-dependence of the SANS 
intensity in CeCoIns can be explained by the strong para- 
magnetic effect, and suggests [i ~ 1.7. When /i ~ 1.7, 
H C 2 is about 38% suppressed by the paramagnetic pair 
breaking from the value of no paramagnetic effect, as 
seen in Fig. [T] The strong paramagnetic contributions 
are also considered as an origin of new FFLO phase and 
first order H C 2 phase transition at higher fields. 

The purpose of this discussion was to demonstrate that 
the paramagnetic effect can change the slope of In 1^1,0 1 2 
as a function of H, and that in the case of strong para- 
magnetic effect |fi,o| 2 does not show exponential decay. 
We note that for further nice fitting to the experimental 
data, there is a room to include the effect by the Fermi 
surface anisotropy or by deformations of the vortex lat- 
tice configuration. 



VI. FIELD DEPENDENCE OF 
MAGNETIZATION 

We discuss the paramagnetic effect on the magneti- 
zation curves. In Fig. [5fa), magnetization curves are 
presented as a function of H for various T at fi = 1.7. 
The magnetization M to tai = B — H includes both the 
diamagnetic and the paramagnetic contributions. It is 
seen that M to t a i(-H) exhibits a sharp rise near H C 2 by the 
paramagnetic pair breaking effect, and that M to tai(^) 
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FIG. 5: (color online) (a) Magnetic field dependence of mag- 
netization Aftotai for (j, = 1.7 at T/T c = 0.1, 0.3, 0.5, 0.7, 0.9 
and 1.0 (normal state) in d-wave pairing, (b) Mtotal as a func- 
tion of T 2 at B — 0.01, 0.02, 0.03, • ■ ■, 0.21. (c) H- dependence 
of factor /3(H) at n = 0.02 and 1.7. 



has concave curvature at higher fields, instead of a con- 
ventional convex curvature. These behaviors are seen in 
experimental data of CeCoIn 5 ,— and previous calculation 
at T = 0.4T C 41 

In Fig. EJb), Mtotal is plotted as a function of T 2 
for various B. We fit these curves as Mtotal (T, H) = 
Mo + \(3(H)T 2 + 0(T 3 ) at low T. The slope /3(H) = 
liniT^o d 2 Mtotal /dT 2 decreases on raising H at lower 
fields. However, at higher fields approaching iJ c2 , the 
slope (3(H) sharply increases. Thus, as shown in Fig. 
[5jc), f3(H) as a function of H exhibits a minimum at 
intermediate H and rapid increase near H C 2 by the para- 
magnetic effect when fi = 1.7. This is contrasted with the 
case of negligible paramagnetic effect (/x = 0.02), where 
13(H) is a decreasing function of H until H c %. 

The behavior of (3(H) is consistent with that of 
'y(H), since there is a relation (3(H) cx dj(H)/dH 
obtained from a thermodynamic Maxwell's relation 
d 2 M total /dT 2 = d(C/T)/dB and B ~ H. 41 In Fig. Q] 
we see that for fi — 1.7 the slope of j(H) is decreas- 
ing function of H at low H, but changes to increasing 
function near H C 2- This behavior correctly reflects the 
ii-dependence of (3(H). The rapid increase near H C 2 is 
clearly seen at lower temperatures, compared with the 
results at higher temperatures^ 



VII. SUMMARY AND DISCUSSIONS 

We studied the vortex states in the presence of strong 
paramagnetic effect by selfconsistent quasiclassical cal- 
culations, which can be used for quantitative estimate 
of the vortex states even far from ff c2 . Calculating the 
spatial structure of the vortex states and local electronic 
states, we clarified the paramagnetic effects on the vortex 
core structure. There, the core radius is enlarged and the 
internal field around the core is further enhanced, due to 
the enhanced paramagnetic moments at the vortex core. 

Qualitatively estimating the ii-dependence of low tem- 
perature specific heat, Knight shift, magnetization, and 
FLL form factor, we showed the relationship between the 
^-dependence behaviors and the strength of the para- 
magnetic effect. The specific heat, Knight shift, and 
magnetization show rapid increase near H C 2, due to the 
paramagnetic pair breaking which is eminent at higher 
fields. The anomalous .ff-dependences of FLL form fac- 
tor of SANS experiment 3 - in CeCoIn5 are also explained 
by the strong paramagnetic effect. This reflects the 
paramagnetic vortex core structure, affecting the inter- 
nal field distribution. These theoretical studies of the 
independences help us to evaluate the strength of the 
paramagnetic effect from the experimental data of the 
.H-dependences in various superconductors. For exam- 
ple, analyses for Sr 2 Ru04, TmNi 2 B 2 C and URu 2 Si2 are 
given elsewhere i 50 ' 51 i 52 

As for CeCoIns, our analysis of FLL form factor in 
SANS experiment suggests /i ~ 1.7. This indicates strong 
paramagnetic effect, so that H c2 is about 38% suppressed 
by the paramagnetic pair breaking from the value of no 
paramagnetic effect. On the other hand, when we com- 
pare the .ff-dependence of the specific heat with the ex- 
perimental dataj22 experimental data j(H) ~ C/T also 
show the concave curve at higher fields, as suggested in 
our calculation. However, at lower fields, j(H) is much 
smaller than that expected by our theoretical calculation. 
This discrepancy indicates that experimental data may 
include other additional contributions, such as antiferro- 
magnetic fluctuations, or ii-dependent bulk properties 
other than the conduction electrons. Therefore, in or- 
der to understand the ^-dependences of the vortex state 
properties in CeCoIns, we need further careful studies by 
the collaboration of experimental and theoretical studies. 
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